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Abstract 

We consider an 0{N) version of a massive, interacting, chiral supersymmetry 
model solved exactly in the large limit. We demonstrate that the system 
approaches a stable attractor at high energy densities, corresponding to a non- 
perturbative state for which the relevant field quanta are massless. The state is 
one of spontaneously broken O(A^), which, due to the influence of supersymme- 
try, does not become restored at high energies. Introducing soft supersymmetry 
breaking to the Lagrangian results in scalar masses at the soft breaking scale rUg 
independent of the mass scale of supersymmetry /i, with even smaller masses for 
the fermions. 



1 



1 Introduction 



Supersymmetry[|T], Q at finite temperature has a number of interesting properties]^, 
H) H) 0' HI- Foremost is tlie appearance of a massless fermion mode, the so-called 
Goldstino. As finite temperature corrections are different for particles obeying 
Bose-Einstein and Fermi-Dirac statistics, a finite temperature state is not invari- 
ant under supersymmetry transformations. Since the transformation parameter 
of supersymmetry is a Grassmann variable, the breaking of supersymmetry there- 
fore implies the existence of a massless Goldstone fermion. 

One situation where such properties might be particularly relevant is during 
the reheating stage after inflation. It has been shown that light fermions may 
play an important role in the reheating process]^, and so it is important to know 
whether massless or nearly massless fermions might be a requirement in super- 
symmetric models due to the phenomenon of the Goldstone fermion. However, 
there is an important distinction between finite temperature physics and the far 
from equilibrium, finite energy density situation relevant to inflation. 

Despite the rich activity in out of equilibrium field theory, the case of su- 
persymmetric models has not yet been properly investigated. In this letter, we 
address a supersymmetric model including all of the superpartner degrees of free- 
dom, including the fermions. The model is one with a global 0{N) symmetry 
which we solve exactly in the limit of large A^0. 

Working at large finite energy density, but far from thermal equilibrium, we 
study the dynamics of the model and see that the system indeed evolves to an 
attractor state for which a set of fermion modes becomes massless. We argue 
that this state is related to the spontaneous breaking of the 0{N) symmetry 
and that, therefore, the fermions are massless because they are superpartners of 
Goldstone bosons. Somewhat surprisingly, the 0{N) symmetry is not restored at 
high energy as would normally be expected [1T0|, As the energy is increased, 
the contributions from the various superpartners to the effective masses of the 
particle modes cancel each other, a result due to supersymmetry. 

We stress that although there are 0{N) symmetric vacua in this theory which 
are degenerate in energy with the 0{N) broken vacua, at high energies the sys- 
tem invariably evolves to the 0{N) broken vacua for which there are massless 
fermions. In this way, the final vacuum state of the system is predetermined 
by the high energy evolution. The possibility that there may be selection rules 
between degenerate vacua resulting from the early evolution of the universe is a 
primary result of this work. 

The model we study contains an 0{N) singlet superfield, characterized by 
the non-zero vacuum expectation value of the scalar component, and an 

0{N) vector supermultiplet. The 0{N) vector fields have time dependent masses 
determined by an order parameter m^(t) reflecting the internal dynamics of the 
system. If (f){t = 0) is far enough from the supersymmetric vacua, i.e., if the 
initial energy density is sufficiently large, the system evolves in such a way that 
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the order parameter (t) asymptotically vanishes, corresponding to zero masses 
for all of the fields of the 0{N) vector multiplet. The state is indicative of the 
spontaneous breaking of the 0{N) symmetry. 

We note that while the ground state upon which this finite energy state is 
built is supersymmetric, the energy is distributed differently among the fermions 
and bosons due to their differing statistics. This is analogous to what is known 
from equilibrium studies at finite temperature P, In a cosmo logical 

context it is expected that the quanta in the highly excited state would become 
diluted with expansion and the universe would eventually find itself approaching 
the underlying ground state. 

It is also worth mentioning that the state is stable to the introduction of small 
soft supersymmetry breaking to the Lagrangian. Such terms explicitly break the 
super symmetry, while not introducing any terms which disrupt the supersymmet- 
ric solution to the hierarchy problem. Here we find that the qualitative behavior 
is unchanged by such terms as the system still evolves to an attractor 0{N) bro- 
ken state. Furthermore, such terms in fact introduce small masses to one set of 
scalar fields of order the soft breaking scale m^, which is taken to be much smaller 
than the overall scale of supersymmetry fi. The fermions meanwhile gain a mass 
several orders of magnitude smaller. 

Although this is only a toy model, the important features of the model - 
the combination of continuous symmetries leading to Goldstone bosons and su- 
persymmetries leading to massless superpartners - are completely general. We 
also note that in models without additional continuous symmetries, such as the 
ordinary Wess-Zumino model, the far from equilibrium system is found to evolve 
toward a state for which the fermions become massless; this case will be stud- 
ied in detail in a future publication [p!^]. The results may be important to our 
understanding of aspects of cosmology, such as supersymmetry based inflation 



models |]13[ and electroweak baryogenesis[Tl, IT6|, IT 



2 The model 



The 0{N) extension of the Wess-Zumino model |T^ consists of a chiral superfield 
multiplet 5*0 = {Aq, Bo]ipo; Fq,Go), which acts as a singlet under 0{N), coupled 
to chiral superfields Si = {Ai, Bf, ipf, Fi, Gi) with i = 1 . . . N , transforming as a 
vector under 0{N). Here, A and F are real scalars, B and G are pseudo-scalars, 
and ip is a Majorana fermion. The superpotential has the form 

WiSo,S.) = ^MS', + -^S', 
^ 1 \ 

We expand in terms of the component fields and eliminate the auxiliary fields 
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via their equations of motion. To allow for a proper large N limit, the expectation 
value of Aq must be of order We therefore set 



(Ao) = VN<P , (So) = . 



(2) 



The field Bq is taken to have zero expectation value for the sake of simplicity. 
We assume that the initial state satisfies the 0{N) symmetry which requires 
that (Ai) = (Bi) = 0. It is convenient to take full advantage of this symmetry 
to define fields A, B, and ip such that J2iAAi = NA'^, EiBiBi = NB'^, and 
J2i '^ii^i = Nipip. The resulting Lagrangian to leading order in is 

— = -d^(j)df'(j)+-d.Ad^A + -d.Bd^'B 
A 



I (/. + X<j>f {A^ + B-) 
'm<P + + 1a [a^ - B^)] [a^ - B 



i - 1 - 1 - 



(3) 



The system may be completely characterized by the equation of motion for 
the mean field 6 



M + K(f) o , , 1 

(p + --^ml + Xm^A^ + B^) + -A(^^) = 

A 2 

and by the time dependent masses of the ip, A, and B fields: 







(4) 



m. 



H + X(j) , 

2 I 2 
2 2 

- m_ 



with. 



m 



X 



M4> + ^k4>^ + ^X{a^-b^) 



(5) 
(6) 
(7) 

(8) 



The expectation values, (ipip), (A"^), and (5^) are determined in the usual way 
in terms of the non-equilibrium mode functions for the individual fields. General 
details, including the renormalization procedure may be found in Refs. 
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Specific details for this model will be provided in future work||12||. 

Note that the appearance of [A^ — B"^) on the right hand side of the expression 
for m?_, Eq. means that this expression plays the role of a gap equation which 
must be satisfied by the dynamics. Also note that the sum rule m\ + m\ — 2m'^ = 
is automatically satisfied for all times as required by supersymmetry. 

Through use of the equations of motion, it is straightforward to show that the 
variation of the Lagrangian (RI) vanishes under supersymmetry transformations up 
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Figure 1: Zero mode evolution; the parameters are /i = 1, M = 4, k = 1, A = 2, 
and 0(0) = 1.5. All masses are scaled by /i which is arbitrary. 

to a total derivative. To this order, the Lagrangian is completely supersymmetric 
with (j) acting ClS db classical background field. 

We stress that the evolution determined by equation @ and the time depen- 
dent masses (^) - (H) exactly solves the quantum field theoretical system described 
by the superpotential (|l|) in the N ^ oo limit to all orders in perturbation the- 
ory. We have made no further approximations and the only simplifications, eg., 
(-Bo) = (Ai) = {Bi) = 0, correspond to choices of initial conditions. 

3 Numerical results 

Finite energy density is imposed via out of equilibrium initial conditions for 
the zero mode 0. The evolution toward the non-perturbative attractor state is 
depicted in Fig. with the corresponding masses of the fields in Fig. |^. We note 
the following characteristics. The evolution begins with large oscillations of 
over the entire classically allowed range of evolution. During this initial period, 
the field fluctuations (A^) and {B"^) grow. After a relatively short period of time, 
the mean field settles down precisely to the point = —fi/X. The result is that 
the N fermions become massless. 

Interestingly, the scalar fiuctuations continue to grow until a state is formed 
for which the fields A and B are massless as well (Fig. §). This configuration 
remains completely stable. We also find that this behavior persists up to energy 
densities much higher than any natural scale in the problem, indicating that there 
is no symmetry restoration and that this represents the generic behavior of the 
system. This is made possible by the cancellation of contributions from (A^) 
and (5^) in the gap equation (|]), a direct consequence of supersymmetry, which 
allows these field fluctuations to become arbitrarily large while not changing their 
contributions to the masses of the field quanta. The growth of these fiuctuations 
also provides the mechanism for driving the fermion mass to zero, see Eq. (^. 
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Figure 2: The effective field masses squared, m\ (top), m% (middle), and 
(bottom) with the same parameters as in Fig. |1]. Each effective mass vanishes at 
late times. 



In order to understand the supersymmetric field configuration reached by the 
non-equilibrium time evolution, it is illuminating to study the effective potential 



for static field configurations as a function (f). It is obtained (see, e.g. |jTO[) by 
maximizing with respect to the potential 

V^(0,m^_) = f M0 + ^ - ^ 



+ 



647r2 



2 " / A 2A 

,2\ I „/_2 



g{m'^) + g{mj,) - 2g{m'^) , (9) 



where g{mf) = mfilnimf / rn?) — 3/2). 

The effective potential V{(j)) is plotted in Fig. ^ along with the tree level 
potential. We see that the effective potential has a minimum at = — yu/A, 
and one finds that at this point all the masses vanish. To understand this new 
minimum, it is helpful to remember that, if one allows for states which break 
0{N), such that, for example, (Ai) = a ^ then one finds 0{N) breaking 



minima at A0 = — /i and Act = ±y2/iM — Kjjfi/X. As a result of the convexity 
theorem, the exact large N effective potential must be flat between these minima 
as in a Maxwell construction. Hence, there must be a new minimum in the full 
effective potential as shown in the figure. 

The state that is reached by the evolution is therefore a phase consisting of 
different spontaneously broken 0{N) states in coexistence. 

The values of and the mass parameters obtained at late times in the out of 
equilibrium evolution correspond precisely to such a state. The state itself is a 
highly excited one and is time dependent in a coherent way. Furthermore, we find 
that the system is attracted towards this configuration once the energy density 
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Figure 3: The large N effective potential (solid line), showing the 0{N) symmet- 
ric ground states at = and = — 2M/ n = —8 and the spontaneously broken 
0{N) coexistence phase at = — /i/A = —1/2. The inset frame magnifies the 
region around the 0{N) broken minimum. In each case, the tree level potential 
(dashed line) is shown for comparison. The parameters are as in Fig. |l|. 

is sufficiently high. 

As a side note, we point out that the effective potential in the direction is 
non-convex. This has to do with the fact that acts as a simple classical zero 
mode whose ffuctuations are completely negligible in the leading order large N 
limit. 

Next, we introduce soft symmetry breaking to the 0{N) model via a scalar 
mass rUs for the A and B fields such that m\ = + + and = 
— + m^. Fig. | shows the result for a value = 10~^. We see 

again that the system reaches the attractor state, and the explicit soft symmetry 
breaking terms provide a mass for the A field equal to \/2ms, while the B field 
remains massless. The fermion mass is found to be three orders of magnitude 
smaller. 

We make the following conclusions. First, supersymmetry may play a very 
important role in the dynamics of the early universe beyond ordinary model 
building. The requirement of massless fermions appearing in the spectrum, in 
particular, may be important to inflationary and reheating dynamics and could 
also play a significant role in baryogenesis. Furthermore, the constraints of su- 
persymmetry can lead to a preferential choice between multiple degenerate vacua 
through the existence of attractor states at finite energy density. Such issues are 
deserving of further study. 
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Figure 4: The effective field masses squared m\{t), m%{t), and m^(t) in the 0{N) 
broken phase including soft masses for the A and B. The parameters are ^ = 10 
TeV, M = 40 TeV, = 100 GeV, « = 1, A = 2, 0(0) = 15 TeV. The late time 
values correspond to uia — V^rUs — 140 GeV, uib — 0, and — 140 MeV. 

References 

[1] D. Bailin and A. Love, Supersymmetric Gauge Field Theory and String 
Theory, (Institute of Physics Publishing, Bristol, 1994). 

[2] M. F. Sohnius, Phys. Report 128, 39 (1985). 

[3] A. Das and M. Kaku, Phys. Rev. D18, 4540 (1978). 

[4] D. Boyanovsky, Phys. Rev. D29, 743 (1984); H. Aoyama and D. Boyanovsky, 
Phys. Rev D 30, 1356 (1984); D. Boyanovsky, Physica 15D, 152 (1984). 

[5] S.-P. Chia, Phys. Rev. D33, 2481 (1986); ibid, 2485 (1986). 

[6] A. Das, Physica A 158, 1 (1989). 

[7] J. P. Derendinger and C. Lucchesi, Nucl. Phys. B 536, 483 (1998). 

[8] D. Boyanovsky, M. D'Attanasio, H. J. de Vega, R. Holman and D.-S. Lee, 
Phys. Rev. D52, 6805 (1995); J. Baacke, K. Heitmann, and C. Patzold, 
Phys. Rev. D 58, 125013 (1998); P. B. Greene and L. Kofman, Phys.Lett. 
B448, 6 (1999). 

[9] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz, P. R. Anderson, Phys. 
Rev. D50, 2848 (1994). 

[10] D. Boyanovsky, H. J. de Vega, R. Holman, and J. Salgado, Phys. Rev. D 59, 
125009 (1999). 



8 



[11] related cases of symmetry non-restoration include S. Weinberg, Phys. Rev. 
D9, 3357 (1974); T. Ross, Phys. Rev. D54, 2944 (1996), G. Bimonte and 
G. Lozano, Nucl. Phys. B460, 155 (1996), J. Orloff, Phys. Let. B403, 309 
(1997), M. B. Gavela, O. Pene, N. Rius, S. Vargas-Gastrillon, Phys. Rev. 
D59, 025008 (1998). 

[12] J. Baacke, D. Cormier, H.J. de Vega and K. Heitmann, in preparation. 

[13] D. H. Lyth and A. Riotto, Phys. Rep. 314, 80 (1999). 

[14] S. Myint, Phys. Lett. B287, 325 (1992). 

[15] G.F. Giudice, Phys. Rev. D45, 3177 (1992). 

[16] J.R. Espinosa, M. Quiros and F. Zwirner, Phys. Lett. B307, 106 (1993). 

[17] M. Laine and K. Rummukainen, Nucl. Phys. B535, 423 (1998). 

[18] J. Wess and B. Zumino, Phys. Lett. B49 , 52 (1974). 

[19] D. Boyanovsky, H.J. de Vega, R. Holman, D.S. Lee and A. Singh, Phys. Rev. 
D51, 4419 (1995). 

[20] J. Baacke, K. Heitmann, and C. Patzold, Phys. Rev. D55, 2320 (1997). 



9 



